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PORTABLi:  GRID  (riB/UV-PORT) 


1 . Introduction 

\ 

The  basic  analysis  capability  is  the  Fields-by-Information  Blending 
(FIB)  methodology  ^developed  by  Meteorology  International  Incorporated. 

<-9  The  FIB  technique  has  been  used  successfully  in  many  atmospheric  and 

oceanic  applications  to  date.  These  applications  include  the  analysis 

f. . ‘Q*C-  3 

of  sea-surface  temperature,  sea-level  pressure,  surface  winds,  and  ocean 
thermal-structure  parameters.  It  has  also  been  used  in  the  diagnosis  of 
clear-column  radiances  estimates  from  side-to-side  scans  of  observed 
satellite  radiances,  for  blending  clear-column  satellite  radiance  information 
with  the  mo-  . conventional  atmospheric  mass-structure  model  parameters 
and  for  oceanic  depth-temperature  profile  retrievals  from  a specified  mixture 
of  upper-ocean  thermal  structure  parameters  and  deep  ocean  climatology. 

In  adcition,  a general  scalar  FIB  capability  has  been  formulated  for  analysis 
of  weather  elements  (e.g.,  fog,  cloudiness)  and  other  scalar  fields.  These 
capabilities  are  documented  in  the  list  of  references  at  the  end  of  this 
report . 

The  general  applicability  of  the  FIB  methodology  results  from  the 
inherent  properties  of  the  scheme — all  information  elements  relevant  to 
the  spatial  distribution  of  an  object  parameter  may  be  brought  to  bear  upon 
the  resultant  analysis  in  proportion  to  their  information  content  relative  to 
the  resolution  and  specific  objective  of  the  analysis.  Hence,  object 
parameter  estimates  and  structure  or  pattern  information,  estimates  of 
parameters  related  to  the  object  parameter  via  linear  relationships  and 
linear  constraints  on  the  resultant  object  parameter  analysis  may  be 
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accommodated  by  the  scheme.  These  estimates  may  come  from  current  and 
near-current  reports,  current  related  analyses,  and  earlier  analyses, 
climatology  or  forecasts. 

In  addition,  the  TIB  methodology  is  conceptually  open-ended  in 
its  ability  to  accept  new  information  elements;  and  it  is  general  in  its 
application  to  grids  of  varying  dimensions  and  resolution  down  to  that 
minimum  mesh  length  afforded  by  the  input  object  parameter  density  and 
distribution . 


-2- 


o 


I'm  mutations 


- — — — The  surface  wind  analysis  capability,  FIB/UV,  has  been  designed 
to  use  a concurrent  sea-level-pressure  analysis  via  a diagnostic  balance 
relationship  for  first-guess  information  for  both  the  winds  and  the  spreading 
of  the  winds.  The  sea-level-pressure  analysis  capability  TIB/P  has  been 
designed  to  provide  this  first-guess  information  for  FIB/UV  in  that  both 
analyses  are  done  on  a staggered  grid, (see  Fig.  l)^.where  the  grid  locations 
for  u , v and  p have  been  chosen  for  optimum  use  of  the  diagnostic 
relationship. ^In^addition , use  of  this  grid  permits  ordering  of  the  unknown 
resultant  object  parameters  so  that  the  line  successive-over-relaxation 
(SOR)  technique  may  be  used  for  solving  the  blending  equations  in  both 
TIB/UV  and  FIB/P. 

A detailed  description  of  the  solution  of  the  blending  system  of 
equations  by  line  SOR  is  included  as  Appendix  A.  Programming  details 
are  given  in  Appendix  B.  Although  expressed  with  respect  to  blending  for 
P* , these  details  apply  equally  to  blending  for  U*  , V*  . 

Consequently,  the  FIB/UV  and  FIB/P  analyses  may  be  run  in  tandem — 
each  FIB/UV  analysis  is  preceded  by  a FIB/P  analysis.  This  increases  the 
information  yield  in  FIB/UV  in  several  respects.  Firstly,  a good  first  guess 
for  the  spreading  parameters  in  TIB/UV  is  produced  by  using  a concurrent 
pressure  analysis  of  equivalent  resolution.  Secondly,  it  permits  the 
assimilation  of  information  available  in  the  pressure  portion  of  the  reports. 
This  information  would  be  difficult  to  incorporate  directly  into  the  wind 
analysis.  Thirdly,  kinematical  extrapolation  in  time,  of  sea-level-pressure 
fields  for  first-guess  information  to  the  sca-lcvel-pressure  analysis,  is 
more  valid  than  kinematical  extrapolation  of  surface  wind  fields  foi  first 
guess  to  the  surface  wind  analysis. 
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1 The  staggered  grid  showing  reference  locations  for  the  pressure  (•), 
the  u component  of  wind  ( x ) and  the  v component  of  wind  ( o ) . 

The  ordering  of  the  grid  points  in  HB/P  is  shown  by  the  solid  sawtooth 
line  and  in  PIB/lJV  by  the  dashed  sawtooth  line.  An  arbitrary  reference 
module  is  outlined  with  dotted  lines.  When  referring  to  locations  on 
the  grid,  the  (Jt.’fl)  grid  notation  is  used.  When  referring  to  the 
ordering  of  the  points  for  blending,  the  (k,n)  notation  is  used. 


Schematics  depicting  the  flow  of  operations  in  FIB/P  and  FIB/UV 
are  shown  in  Figures  2 and  3,  respectively.  The  information  elements  for 
FIB/P  and  TIB/UV  are  shown  in  Figures  4 through  7.  Details  of  the  FIB/P 
and  FIB/UV  analysis  capabilities  are  given  in  references  [31]  and  [27]. 
The  present  discussion  will  be  limited  to  those  adaptations  necessary 
for  application  on  the  portable  grid. 
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Fig.  3 Schematic  Flow  of  Operations  for  FIB/UV 
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Fig.  4 The  arbitrary  FIB/P  l ,n  module  for  the  staggered  grid  consisting 
of  one  corner  point  one  center  point  ( iH/2  ,m  +1/2  ) , two 

sides  and  the  interior  area.  Reference  points  within  the  module 
are  shown  for  the  pressure,  p , and  associated  weight,  A . 


Fig.  5 Reference  locations  shown  for  the  TIB/P  first-il ifference  spreading 
parameters  defined  below.  The  associated  weights  are  given  by 
the  adjacent  upper  case  letters.  The  indices  l and  ^ assume 
both  integer  and  half-integer  values. 
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Reference  locations  shown  for  the  TIB/P  second-difference , 
Laplacian  and  the  cross-difference  spreading  parameters  defined 
below.  The  associated  weights  are  given  by  the  adjacent  upper 
case  letters.  The  indices  l and  n assume  both  integer  and 
half-integer  values  . 
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Fig.  7 The  arbitrary  l,"i  area  module  (consisting  of  one  corner  point,  two  sides,  and  the 
Interior  area)  and  adlolnlng  modules.  Reference  locations  In  the  module  are  shown 
for  the  u and  v components  of  the  wind  velocity,  and  the  finite  difference  parameters: 
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The  weights  of  the  values  are  given  by  the  adjacent  upper  case  letters  In  the  figure. 
The  reference  location  for  the  pressure,  p , from  TIB/P  are  also  shown. 


(denorali/.ation  of  the  Grid  Location,  Orientation,  Sl/.o  and  Mesh  Length 


The  analyses  are  done  on  a polar  storcocjraphic , rectangular  grid. 
The  grid  dimensions  are  variable  up  to  (MxN)  = (63 xb3)*  and  the  grid 
mesh  length  is  variable  down  to  the  resolution  afforded  by  the  density 
and  distribution  of  wind  and  pressure  reports  in  the  area  of  application. 

A complete  specification  of  the  portable  grid  includes  specification 
of  the  center  of  the  grid,  the  rotation  relative  to  the  hemispheric  grid,  the 
dimensions  M and  N , and  the  grid  length  as  a fraction  of  the  standard 
63x63  hemispheric  grid  length  (i.e.,  the  grid  length  factor).  Figure  8 
shows  the  standard  Nil  63x63  grid  with  the  embedded  greater  Mediterranean 
regional  grid.  Table  1 specifies  the  grid  parameters  for  each.  The  standard 
TNWC  63x63  hemispheric  grid  has  a mesh  length  of  381  km  at  60°  N.  The 
regional  grids  used  for  FIB/SLP-MED  [5]  and  FIB/UV-MED  [13]  employ  a 
grid  length  factor  of  0.25. 


The  total  number  of  grid  points  on  the  MxN  grid  is  MxNx2  including  the 
module  center  grid  points.  See  Tig.  1. 
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Fig.  8 TNWC  Northern  Hemisphere,  polar  stereographic  63x63  grid 
with  embedded  greater  Mediterranean  region  indicated. 
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Table  1 


Grid  Parameters  for  the  FNWC  Northern  Hemisphere, 
Polar  Stereographic  63x63  Grid  and  the  Embedded 
Greater  Mediterranean  Grid 


NH 

MED 

Center  of  Grid: 

1 = 31.  1 = 31 

1 = 46,  J = 33  (on  NH  Grid) 

Counter-Clockwise 

Rotation: 

0° 

0° 

MxN: 

63x63 

63x63 

Grid  Length  Factor: 

1.0 

.25 
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Generation  of  Term  in  Parameters 


Although  the  areas  of  primary  interest  arc  over  the  oceans,  it  is 
anticipated  that  the  specific  grid  encompassing  many  of  these  regions  will 
include  adjacent  land  masses  and  islands.  Consequently,  the  terrain 
parameters  used  in  the  FIB/UV  and  FIB/P  analyses  will  bo  required  for  each 
new  portable  grid  employed  which  includes  land  masses  and  islands. 

These  parameters  are  a function  of  the  land  topography  and  the  relative 
percentages  of  land/sea  within  the  data  assembly  and  reference  modules 
on  the  grid. 

In  the  FIB/UV  analysis,  the  weight  associated  with  the  assembled 
wind  report  value  in  each  data  assembly  module  is  reduced  as  a function 
of  the  module  terrain  roughness  and  land/sea  mix--both  factors  contributing 
to  non-representativeness  in  the  winds  as  opposed  to  modules  over  the 
open  ocean  or  flat  land.  In  addition,  over  rough  terrain  the  winds  between 
adjacent  grid  points  will  be  less  strongly  correlated  than  over  smooth 
terrain.  Therefore,  first-difference  and  other  higher-order  spreading  weights 
are  reduced  as  a function  of  the  specific  reference  module  terrain  roughness. 
These  formulations  arc  given  in  more  detail  in  reference  127]. 

In  the  riB/P  analysis,  the  weights  associated  with  the  assembled 
pressure  report  values  are  reduced  over  land  as  a function  of  mean  assembly 
module  elevation  to  reflect  the  uncertainty  in  the  reduction  of  pressure  to 
sea  level.  In  addition,  assembled  first-difference  weights,  which  reflect 
the  assimilation  of  wind  reports  via  the  balance  relationship,  are  reduced 
as  a function  of  mean  assembly  modulo  elevation  and  terrain  roughness. 

The  first  reduction  accounts  for  extrapolation  to  sea  level  and  the  second 
for  decreased  correlation  between  points  because  of  uneven  terrain.  Also, 
the  weight:;  ol  higliei  -oidei  .pleading  paiamelei:;  aic  i educed  because  ol 
uneven  terrain.  These  [ormulnlions  are  given  in  more  detail  in  reference  [31]. 
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The;  len.iin  parameters  arc  presently  derived  f rom  the  I'NWG, 
Northern  Hemisphere,  hiyh-resolution  Lt ;r i <i in  dale  available  <ji  1 magnetic 
tape.  The  capability  to  process  this  raw  data  into  the  required  terrain 
parameters  for  arbitrary  reference  ami  data  assembly  modules  has  been 
developed  under  earlier  contract  with  FPHF  [30].  The  present  capability 
has  been  extended  to  compute  all  the  required  terrain  parameters  necessary 
for  both  FIB/UV  and  TIB/P  on  the  specific  grid,  to  output  and  catalogue 
(for  future  reference)  the  parameters  in  the  required  format  to  interface 
with  FIB/UV  and  TIB/P,  and  to  generalize  the  formulation  as  a function 
of  grid-mesh  length  for  assimilating  the  parameters  into  the  analyses. 

Grids  extending  into  the  Southern  Hemisphere  will  not  have  the 
benefit  of  the  terrain  parameters  since  the  basic  terrain  data  for  the 
Southern  Hemisphere  are  not  currently  available.  However,  the  capability 
to  use  the  terrain  parameters  will  be  in  the  models  whenever  the  Southern 
Hemisphere  data  become  available. 


Generalization  of  Rccvaluation  Procedures 


In  the  TIB  analysis  scheme,  three  cycles  are  employed.  After  the 
first  two  cycles,  input  data  weights  are  reevaluated  as  a function  of  the 
consistency  between  individual  input  values  and  all  other  information. 
The  weights  of  individual  information  elements  which  arc  not  consistent 
enter  the  second  and  third  data  assembly  cycles  with  reduced  weight 
according  to  the  degree  of  disparity.  The  formulation  of  this  reevaluation 
is  dependent  on  the  resolution  of  the  specific  grid,  as  is  the  initial 
weighting  of  ihe  Information  elements. 
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( '.(Minal  ix.aliun  of  Weighting  Formulations 

liacb  clement  of  input  information  in  a FIB  analysis  is  weighted 
according  to  its  reliability  (see  Figs.  2 through  7 for  the  specific 
parameters  and  associated  weights  used  in  F1B/P  and  FIB/UV).  In  general, 
the  weights  are  a function  of  the  grid-mesh  length  (i.e.  , resolution)  as 
well  as  of  the  quality  of  the  information.  In  addition,  the  assembled  weights 
of  first-difference  pi  assure  values  in  riB/P  are  reduced  before  blending  to 
reflect  the  inherent  deficiencies  in  the  diagnostic  balance  relationship 
used  to  express  reported  winds  as  equivalent  pressure  differences.  This 
weight  reduction  is,  among  other  things,  grid-mesh-length  dependent. 


First-Guess  Field  Preparation 

The  first-guess  field  for  FIB/UV  is  prepared  directly  from  the 
concurrent  TIB/P  analysis  using  the  diagnostic  balance  approximation. 

The  first-guess  pressure  field  at  analysis  time  is,  in  general,  a 
combination  of  the  kinematically  extrapolated  fields  from  the  previous 
analysis  and  a later  analysis.  When  only  the  previous  analysis  is 
available,  no  combining  is  necessary.  Kinematical  extrapolation  is 
performed  with  geostrophic  winds  derived  from  the  SR500  field  height 
values . * 

This  method  of  computing  the  first-guess  field  provides  for  both 
forward  and  backward  time  continuity  in  the  TIB/P  analysis  model,  and, 
consequently,  in  the  concurrent  FIB/UV  analysis. 


* 

The  SIGOO  field  ir;  the  actual  500-mb  height  field  with  the  disturbance 
scale  removed. 
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Sol  -up,  St • irt  — u|>  and  Schodul  ing  Options 

The  tandem  riB/UV-FIB/P  analysis  capability  available  on  an 
"as-roquired"  basis  on  a portable  qrid  and  with  provision  for  backward 
and  forward  continuity  in  time  is  a powerful  and  flexible  tool.  This 
flexibility  requires  the  consideration  and  specification  of  several  types 
of  options . 

1 --  Set-up  options:  Before  the  model  is  run,  certain  tasks  must 

be  performed.  These  include  specification  of  the  grid 
parameters  and  generation  of  the  required  terrain  parameters 
as  discussed  above. 

2 --  Start-up  options:  The  model  must  be  run  for  some  specified 

"start-up"  period  preceding  the  "operational-use"  period. 

3 --  Scheduling  options:  The  actual  run  times  for  the  analyses 

must  be  specified. 

Items  2 and  3 above  will  now  be  discussed  in  more  detail. 

Start-up  Options 

Each  new  request  to  begin  a series  of  analyses  will  have  to  specify 
an  operational  use  period  for  the  analyses.  The  end  points  of  the 
operational  use  period  will  be  given  as  t,^  and  t^  . The  operational  use 
interval,  then,  is  ^ ~ l2  > 0*  The  start-up  time,  tQ  , for  the  model  will 
have  to  be  sufficiently  in  advance  of  t2  so  that  the  effect  of  the  initial 
conditions  introduced  by  the  bootstrap  technique  have  been  eliminated. 

We  will  indicate  this  interval  for  elimination  of  bootstrap  effects  by 
t|  - t(|  whore  tj  ' t . In  addition,  I ho  ana  lyse:;  wi  1 1 bo  run  at  fixed 
inleivjls,  At.  those  telal  ion:. hips  ate  summ.u  i/.od  in  fig.  'J . 
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Fig.  9 Schematic  of  Time  Intervals  and  Relations. 


Due  to  the  provision  for  forward  and  backward  continuity  in  FIB/P, 
the  start-up  period  should  be  short.  Also,  the  existence  of  a concurrent 
hemispheric  or  regional  sea-level-pressure  analysis  into  which  the 
portable  grid  can  be  embedded  will  further  shorten  the  start-up  period. 
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Scheduling  Considerations 


Scheduling  the  runs  into  the  operational  job  stream  involves  many 
considerations.  Among  these  are  the  following: 

1.  Operational: 

a.  Requirement  to  interface  with  other  products  (c.g.,  input 
to  forecast  or  analysis  model). 

b.  Requirements  to  transmit  to  the  field. 

2 . Da  ta : 

a.  Lag  time  between  observations  and  receipt  of  observations. 

b.  Data-rich  vs.  data-poor  areas  and  periods. 

c.  Degree  of  updating  of  analyses  to  be  performed. 

3.  Time  Resolution: 

a.  Analysis  period. 

b.  Extent  of  forward  and  backward  continuity  desired. 

4.  Processing  time: 

a.  Amount  of  computer  time  available. 

b.  Period  when  computer  time  available. 

c.  Model  running  time. 

5.  Historical  study  vs.  operational  use. 

The  scheduling  for  a specific  application  will,  in  general,  involve 
all  of  the  above  considerations.  No  single  scheduling  will  be  optimum  to 
moot  all  possible  applications  of  the  model. 
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3. 

; iupplementi  > 1 Details 

1 JB/P 

3.  1 

QB/ 1 1 < ■'  mipi  >nonl  < )pcr 

utinrir; 

1.  Sot  Constants:  analysis  time,  grid  definition  and  identification,  and 
analysis  cycle  time  increment. 

2.  Read  elevation  Data. 

Set  to  zero  if  not  found  or  not  yet  produced. 

3.  Set  Adjustable  Constants. 

4.  Read  First-Guess  Analyses. 

a.  Latest  SR500  analysis  for  kinematic  extrapolation 

b.  Previous  and  following  FIB/P  125x125  analyses  j previous  time=IAII(l)  •'  0 

nearest  to  current  time  (up  to  Jd)  hrs)  1 following  time=TAH(2)  > 0 

Use  riB/SLP  version  if  and  only  if  no  TIB/P  version  analyses  found. 

5.  Determine  Number  of  Relaxation  Passes. 

N TOT  HR  = 1AH  (2)  - IAH(l)  I IAH(m)  | 5 9 

a.  If  only  previous  or  following  analyses  found 
NPPY  = NTOTHR/2  4-  1 

b.  If  both  found 

NPfY  = NTOTHR/3  - 1 


Aces  of  First  Guesses 

N PTY 

Number  of  Passes 

1 found 

2 found 

Cycle  1 

Cycle  2 

Cycle  3 

6 

6,9 

5 4 

15 

10 

6 

6,6 

3 

13 

9 

6 

3 

3,6 

2 

10 

8 

5 

JL1 

1 

8 

7 

5 
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(>.  Gombino  will)  Tropical  l’iold(s). 

W 2 .-.in  0 (II  < W < \ ) O latitude 

^ W - weight 

II  hemispheric  (FIB) 

T II  T = tropical  (global  band) 

INTRP  used  in  tropical  field. 

7.  Kinematically  Fxtrapolate  Combined  Field (s)  in  3 hour  steps. 

FKX  = 0.5B  x SR500  geostrophic  velocity 

a.  Linear  interpolation  used  in  moving  l,j  indices  with  SR500 . 

b.  STGNTR(12  point)  interpolation  uses  these  I,  I \alues  to 
interpolate  in  previous  and/or  following  analyses. 

8.  Combine  Fxtrapolatcc!  Fields  into  First  Guess,  inversely  weighted 
by  time  (IAI1) . 

9.  Read  Observational  Data. 

Location,  age,  type,  name,  wind,  pressure  and  pressure  tendency, 
ship  movement , classification,  etc. 

a.  ships  (hourly) 

b.  synops  (3  hourly) 

c.  airways  (hourly) 

10.  Sort  Land  Reports  north  to  south. 

Save  nearest  to  0 hr  when  duplicates  found. 

11  . Move  T3  hr  reports  back  to  0 hr  using  tendency  group  (ship  and  synop) 
and  movement  code  (ship). 

12.  Sort  Ships  by  Name. 

Flelete  same  name  report  if  within  distance  eguivalent  to  2 latitude 
and  pi  c run  < • williiu  Xmb  ot  age  / 0. 
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13.  Read  Artificial  ("Bogus")  Data, 
a . I'oiccil  pnssuie/wind  icporl 

b.  Forced  accept  or  reject  of  wind  or  pressure  by  station  name  or  area 

c.  Tropical  storm 


H.  Sort  all  Data  by  J (required  for  ordered  assembly). 

15.  Form  Final  Data  List  (4  words/report). 

Determine  report  weights. 

Convert  wind  to  equivalent  pressure  difference. 

16.  Set  Accept/Reject  Bits  in  data  list  from  area  and  name  type  boguses. 

17.  Set  up  Tropical  Pressure  Fields  from  TRO  Boquses — used  for  pressure 
differences  only. 


18.  a.  Compute  Spreading  Parameters  from  first-auess  field 


first  deriv  weights  = Wl 
higher  deriv  weights  = W2 
Laplacian  weights  = W2T 


- 0 where  undefined 


) 


wi  = o.onoi 

W2  =0.02 
/ W2T  = 0.1 


b.  Use  "tropical  pressure  field”  only  in  vicinity  of  tropical  storms 
(vicinity  defined  as  2 x 30  kt  radius). 


19.  Get  Mean-Square  Laplacian  Field  (symmetrical  5-point  mean).  Set 
low  order  bit  only  in  vicinity  of  tropical  storms. 

20.  Add  Terrain  Roughness  and  Laplacian  Variance  to  higher  derivative 
fields,  except  in  vicinity  of  tropical  storms. 

added  variance  = (CS)  o2,  MCL)L2  ( S 1,0 

elev  CL  = 0.6 

2 1 . Get  h igher  der  iv  ( R , S ,T  , W ) portions  of  a , b , c , O'  , ft  , "Y  , 6 , € , R fields 
and  Save. 
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Main  Loop  Bog  ms 

22.  Assemble  Data. 

a.  pressure 

b.  6 first  dcriv  components 

23.  Add  elevation  Variance  to  Pressure  Assembly. 

Added  Variance  = CLP  x (LLLV)^  CLP  = 0.001  x (ft/100)  ^ 

24.  Add  Llovation,  Terrain  Roughness  and  Balance  Approximation  Variance 
to  first-difference  assemblies. 

25.  Compute  a , b,  c , a , B , V , 6 , e , R fields,  omitting  R,S,T,W  terms. 

26.  Combine  Tields  with  and  without  R,S,T,W  terms. 

27.  Initialize  Fields  for  Forward  Elimination,  and  Backward  Substitution. 

28.  Compute  p*  Iteratively  by  line  SOR.  Jump  to  30  if  final  cycle. 

29.  Reevaluate  Pressure  and  Wind  Reports,  then  jump  to  22. 

Skip  this  step  if  final  cycle. 

30.  Output  Analysis,  Varian  Data  List  and  Reject  List. 

Analyses  output  in  125x125  form  (points  added  from  staggered  grid) 
and  63x63  form  ("center"  points  omitted). 
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1.2 


I'll!  P_(  lomponenl  I )( > 1 . 1 Ms 


& 


(2) 

(3) 

@ 

(D 


© 


I I IV  P_KJjiem.il  ii'  lIxli  il'oKijion  (Relet  lo  .Section  i.l  item  7) 

Zoom  steering  field  to  same  scale  as  object  field  with  "center 
points"  omitted.  Not  required  if  doing  hemispheric  analysis. 

Pack  1,J  starting  values  into  steering  field. 

extrapolate  I,J  displacement  fields  in  3-hr  steps.  If  hemispheric 
analysis,  displacement  is  damped  equatorward  from  20°  lat. 

Fill  in  I,[  fields  at  "center  points"  using  4-point  means. 

Interpolate  at  I,J's  in  starting  field  for  all  grid  points. 

If  regional  analysis,  interpolate  in  hemispheric  field  when  I,J  is 
closer  than  "x"  grid  lengths  from  edge.  Combine  the  values 
obtained,  weighted  by  distance  from  edge: 


W = , W = 1 - W , W =wu  V W u. 

hem  x obi  hem  comb  hem  obj 


1 


© 


where  0 
x 


W. 

hem 


1, 


1.5 

GLf 


0.5  (truncated),  and 


d = fractional  number  of  grid  lengths  from  nearest  edge 
(measured  in  direction  normal  to  edge). 


Repeat  (T)  through  (5)  if  extrapolating  backward  as  well  as  forward 
in  time.  Then  combine  these,  weighted  by  amount  each  is  off-time. 


Note,  the  above  method  does  not  require  adding  a border  region  around  the 
regional  grid,  as  previously  done. 
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riH/l>  Pressure  Assembly  (Refer  to  Section  1.  1 item  22a) 


Initial  report  weight  - An  = 0.5 

2 

(tuned  for  approximately  67%  X < 1 ) 
in  each  assembly  cycle: 


p = p + p - p. 

Fn  o rpt  Hi 


where 

P = assembly  ("extrapolated")  pressure  at  ref.  grid  point 
n 

p = first-guess  pressure  at  ref.  grid  point 

0 

p = reported  pressure  at  station  location 
rpt 

p = first-guess  pressure  interpolated  at  station  location 

1 

STGNTR  12-point  interpolation  used 
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I'lll/T  I’it s L- 1 inference  Assembl i o : ; (Refer  to  Section  3. 1 item  22l>) 


r*H 


Definition: 

GLF  = grid  length  factor 


GLSOO 

VWCON 


( OLE ) 

\1  t-  sin  60 °J 


= 0.47*  (GLS60)2 


( - 1 (N.  Hem) 

\ = 0.25  (Mod.) 


Computed  for  Each  Wind  Observation: 


wind  variance  (wind  units): 

wl:l  = 9 f (rr-20)2/40 
FF  = observed  speed  (m/sec) 


report  variance  in  pressure  diff  units: 

VWND  = (VWCON)  [ ( 1 + sin  0 ) ( sin  0 ) ] 2 (WEL) 

report  weight  for  assembly 

Bn  = min-  of  » 7-875 

Six  components  (b,c,d,e,f,g)  are  assembled  with  weight  Bn  at 
grid  points  which  are  dependent  on  exact  location  within  module 
(see  p.  27,  M-209  Technical  Summary). 


The  value  0.47  was  derived  for  c,e  components  but  is  also  used  in 
weighting  the  b,d,f,g  components. 
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FIB/P  I*  i » ; '•  t — I )il  lorotico  Variance  Additions  hy  (aid  | •,  >i  pi :: 
(Refer  l<>  Sectii >n  3 . I i l<  an  2d) 

Variance  added  hy  grid  point  after  assembly  for: 

1.  Elevation 

2.  Terrain  Roughness 

3.  Ba lance  Approximation 

1.  Elevation  Variance  = (CE)  (E) 

E = mean  modal  elevation  in  module  area 
CE  = 0.2  for  E in  ft/100 

2 

2.  Terrain  Roughness  Variance  = (CS)  (cr  ) 

o = std  dev  of  modal  elevations  in  module 
CS  = 1.0  for  a2  in  ft/100 

L 


3.  Balance  Approx.  Variance: 


I'tB/P  Pressure  Reevaluation  (Refer  to  Section  3.1  item  30a) 


X 


2 

P 


2 2 

a + a 

c q 


where 


- report  pressure,  extrapolated 
p*  = analyzed  pressure  at  grid  point 


a = class  variance 

c 

o 

cr  = gradient  variance 

g 

= mean  diagonal  first-difference  squared  / F 

g 

Fg  = GSQF  (GLF)2  GSQT  = 6 

2 

ct  limit  = 4 

g 


Modification  for  low  pressure: 


X 


2 

P 


is  multiplied  by: 


F 

P 


p*  - 920 
100 


(mb) 


0.4  5 rp  5 i.o] 


Reevaluated  weight: 

ANR 


AN/('  '(XP-  >)2) 


ANR  - reevaluated  weight 
AN  = class  weight 
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IWP  I'irst- Difference  Heev< i lug t ion  (Refer  L<  > Section  3 . I item 


whore 


1/B 

n 


Bn  = first  cycle  assembly  weight  (b,c,d,c,f,q  components) 

c,c  = J,I  direction  first  differences 
2 2 

A , = A for  first  difference 

d 

2 

Modified  A for  high  wind  speed: 

2 

X ^ is  multiplied  by: 


3.5 

MIN 

[' 

t *2  *2\l/2 

[u*  + v*  ) , \ 

(2  2 s 

[u  + V ) 

1 r n > 

i1/2] 

Second  Test: 

xd<2’ 


where 


(c  + c*)2  + (e  + e*  )2  + XLWGK 
v n ' ' n ' 


XLMTW  = 6 


TXLWG  = XLMTW/0.7 
XLWGK  = 0 . A (GLF)2 


Rejecl 


if 


X 2 or  X2  (2)  > XLMTW  . 
d d 


30b) 


TXLWG 


-30- 


For  reports  kept: 


reduce  weight 


if  X 


> 


1 : 


B 

r 


B 

n 


1 + (X  - l)2 
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A . Supplemental  Details  FIB/UV 

A . 1 FIB/llV  ( )<  > 1 1 1 p<  mi  Mil  C )pt;i  .1 L ion:; 

Define  portable  grid  information  and  adjustable  constants 
Center  I,J  (in  63x63  grid) 

Scale  factor  (fraction  of  63x63  grid  length) 

Grid  dimensions 

2.  Read  elevation  data 

If  data  not  available,  zero  values  are  assumed 
(correct  over  ocean  areas). 

3.  Read  current  pressure  field  and  diagnose  U,V. 

A.  Compute  K,R  weighted  by  fraction  of  land/sea  coverage  upstream. 

5.  Add  terrain  roughness  variance  to  first-guess  variance  fields. 

6.  Get  first-guess  divergence  and  vorticity  from  U0,VQ  . 

7.  Get  R and  S terms  from  first-guess  fields,  omitting  weights  to 
be  assembled. 

8.  Read  data  and  convert  to  internal  format . Eliminate  duplicates  . 

a.  K , R computed  explicitly  from  observed  wind  and  pressure 
gradient  for  each  observation. 

b.  Special  station"  list  search  for  matches.  Observations  in 
spec.al  list  are  flagged  and  given  listed  weight. 

8.  Main  loop  begins — assemble  data.  First  pass  assembly  weights 
reduced  on  basis  of  first-guess  U,V  components,  when  large 
discrepancies  are  computed,  band  rcfKJils  assembled  first,  then 
ships.  Elevation  variances  added  to  land  assembly  only. 
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10.  Arid  assembled  weights  to  R,S  terms. 

11.  Set  up  matrices  for  line  SOR. 

U.  Get  blended  U,V  fields  by  line  relaxation.  Go  to  #15  if  first 
pass  through  main  loop. 

13.  U.V  values  at  staggered  grid  points  converted  to  standard  grid 
points  by  cubic  interpolation. 

M.  Output  U,V  and  direction,  speed  fields.  Go  to  #16. 

15.  Reevaluate  U.V  components.  List  ship  reports  and  any  stations 

in  special  list.  List  includes  name,  location,  nearest  grid  point, 
reported  and  analyzed  direction  and  speed,  vector  difference,  K 
and  R computed  from  observed  wind  and  pressure  gradient, 

o 

beginning  and  reevaluated  weights,  X and  observed  and  nearest 
grid  point  U and  V components.  A special  notation  is  printed 
for  special  stations  and  for  rejects.  Go  to  #9  if  first  cycle. 

2 

6.  List  X summaries — frequency  distribution  of  X2  values  by 

report  type,  and  mean  and  median  X and  reported  wind  speed 
by  classes  of  analyzed  wind  speed  for  ship  reports. 


4.2  Fib/UV  Adjustable  Constants 


Following  constants  affect  the  grid  placement  and  are  defined  in 
the  program  call: 

1,1  (integer)  define  the  center  point  relative  to  the  standard  grid 
N (integer ) defines  the  grid  dimension 

G (real)  defines  the  mesh  length  as  a fraction  of  the  standard  gi  id 
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The  above  values  are  order  independent,  Example  of  call: 


riBUV  (6 , 3 ,G  0.25,1  -25,  J 41,N=29) 


This  first  parameter  (6)  defines  the  observation  time  as  6 hours  after  the 
current  time  in  the  computer.  The  second  parameter  (3)  indicates  a 
portable  grid  application.  Other  numbers  are  used  for  fixed  areas. 
G=0.25  indicates  a 1/A  mesh  grid.  1=25,  J=4  1 indicates  the  grid  is 
centered  at  this  grid  point  in  the  63x63  grid  numbered  from  0 to  62. 

N=29  indicates  the  dimensions  of  the  grid  are  29x29. 


The  following  constants  arc  defined  in  the  program  under  the 
heading  "Adjustable  Constants" 


XKPL 

XKPS 

RLAND 

RSEA 

VLW 

VSW 

AWT 

DWT 
QWT 
EFWT 
GHWT 
ELS  CL 

SCLRP 

VLSMAX 


land  value  of  frictional  constant  K 
sea  value  of  frictional  constant  K 
land  value  of  frictional  constant  R 
sea  value  of  frictional  constant  R 
class  variance  of  land  wind  report 
class  variance  of  sea  wind  report 
first-guess  A weight  (before  addition  of 
other  variances) 
first-guess  D weight 
first-guess  Q weight 
first-guess  E and  F weight 
first-guess  G and  H weight 
first-  and  second-difference  elevation 
roughness  scaling  faolot 

U ,V  assembly  elevation  roughness  scaling  factor 
maximum  variance  addition  to  11  ,V  assembly 
due  to  land/sea  effect 
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5 . Running  Procedure 

Tho  following  control  cards  arc  required  to  run  a portable  grid 
application.  The  example  shown  is  for  the  current  time  on  a 1/4  mesh 
19x19  grid  centered  at  [31,13]: 

[OBCARD  ,F,C300  ,T300 . 

APLIB(,*FIBP,  *FIBUV) 

RFE.300. 

FIBP(0 , 3 , N=19 , 1=3 1 , J=13  , F=0 .25) 
riBUV  (0,3,  N=19 , 1=3  1 , J=13  , F-0 .25) 

The  first  parameter  in  the  call  indicates  the  number  of  hours 
before  ( - ) or  after  ( + ) the  current  machine  data-time  group.  The  second 
parameter  is  always  3 to  indicate  a portable  grid  application.  The  other 
parameters  are  respectively,  the  grid  dimensions,  the  I and  J location  of 
the  center  in  the  standard  grid  system,  and  the  mesh  size  fraction. 

The  above  application  requires  approximately  19  sec  central 
processor  time  on  tho  CDC  6500  for  FIB/P  and  15  sec  for  TIB/UV. 
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AITKNDIX  A 


Solution  of  tho  Blending  System  of  liquations  by  Line  SOK 


as 


The  system  of  blending  equations  to  be  solved  may  be  written 


CP*  = R 

where  the  ordering  of  the  elements  within  the  vectors  P*  and  R will 

be  that  of  increasing  values  of  k within  increasing  values  of  n (i.e., 

the  k,n  grid-point  notation  introduced  in  Section  1.3  and  shown 

schematically  in  Fig.  3).  This  permits  us  to  partition  P*  and  R into 

subvectors  P*  and  R , respectively,  where  each  subvector  corresponds 
-n  — n 

to  a row  of  unknown  pressure  elements  or  forcing  functions.  These 
vectors  are  summarized  as  follows: 


-Al- 


r 


Referring  to  the  13-point  stencil  for  the  blending  system  of 

equations  (Fig . /),  it  can  be  seen  that  when  the  center  point  (i.'n)  of 

the  stencil  is  the  corner  point  of  a grid  module  (integer  values  of  t.  and 

”2  ),  the  subvector  P*  , , will  have  three  non-zero  coefficients  while 
~n*  1 

subvectors  P*  and  P*  , will  each  have  five  non-zero  coefficients. 

■ n ~ n- 1 

However,  when  the  center  point  ( I ,”  ) of  the  stencil  is  the  center  point 

of  a grid  module  (half-integer  values  of  l and  r),  the  subvector  P* 

~n-  I 

will  have  three  non-zero  coefficients  while  subvectors  P*  and  P*  . will 

-n  ~nH 

each  have  five  non-zero  coefficients . 


-A3- 
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An  equivalent  partitioned  form  of  Cq.  (5)  is  then 


~ — 

“ - 

“ — 

A,  B , 

p* 

R i 

fa  1 1 

-1 

- 1 

B I A „ B „ ( 

ie 

P* 

R _ 

« 1 ^2  £2 

l2 

2 

T 

B A B 

p* 

R 

fan-1  tan  «n 

~n 

— n 

T 

B A , B 

P* 

R ..  , 

, ~ M-2  „M-1  5M-1 

-M-l 

'M-l 

V 

x y t 

B A 

P* 

R 

Sm-i  Sm 

~ M 

where  C has  been  rewritten  as  a block  tridiagonal  matrix.  The  elements  of 

An  come  from  the  middle  row  of  the  stencil  » the  elements  of  Bn  come  from 
~ T 

the  top  row  of  the  stencil  - and  the  elements  of  B,  . come  from  the  bottom 

n-  i 

row  of  the  stencil  . Both  the  A „ and  B submatrices  are  KxK  and 

« n 

pcntadiagonal . In  addition,  the  A are  symmetric  and  positive-definite. 

ts  n 

The  matrix  C may  be  written  as 

C = A + B (A  ] ) 

f ; fj  » 

where 


A 


ding 


A - , A , 

fa  1 M i 


£ M 


(A2) 
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and  B is  given  by 


O 

S3 


B 


1 


B 

£3 


T 

1 


O 


O 


The  blending  system  equations  may  then  be  rewritten  as 


(a  + B j P7 


(A3) 


or  as 


A P*  + B . P*  , + B P * = R 

« n ~n  fan-1  —n-1  » n ~n  + l — i 


n = 1,2, ...  ,M  (A4) 


whei  e 


B 

O 


B 


M 


O 
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equation  (A  l)  may  lx-  solved  lot  1’  as 


* -1 

= A 


P 

a 


R - R , P , 
n ~n  - 1 — n - 1 


R P 

c-l)  II  » 


(A  r> ) 


The  solution  by  line  relaxation  is  defined  by 


(p*  )<rtl)  = 

(p*  )(r) 

+ 

'p*  (r+1/2) 

/p*  i(r) 

UU 

V~n  / 

V — n ) 

- n 

\ n / 

- 

(A  6) 


where  the  superscript  (r)  refers  to  pass  numbe  r and  where 


(£*„  )<r+1/21 


(r) 

(r+1) 


. (A 7) 


The  upper  superscript  (r)  applies  to  the  teim  in  brackets  for  the  first  half 
of  each  pass  in  which  either  the  even  or  odd  numbered  n (hereafter  referred 
to  as  red  lines)  arc  solved  for.  The  lower  superscript  (r+1)  applies  to  the 
term  in  brackets  for  the  second  half  of  each  pass  in  which  the  remaining 
(either  odd  or  even)  n (hereafter  referred  to  as  black  lines)  are  solved  for. 
The  method  converges  for  0 < id  < 2 where  1 < uu  < 2 implies  overrelaxation. 
The  inversion  specified  by  Eq.  (A7)  is  obtained  explicitly  by  a triangular 
decomposition  of  the  A n matrices  into  the  product  of  an  upper  and  a lower 
triangular  matrix,  followed  by  Gaussian  elimination. 

We  will  n^w  present  the  details  of  the  method  of  solution  by  line 
relaxation . 

If  we  let 


*n 


(r  + 1/2) 


(AB) 


-Ah- 


(A1)) 


and 


r 

— n 


(r) 

(r+1) 


(r)l 

(r+D 


then  Eq.  (A7)  may  bo  written  in  the  equivalent  form 


A 

^ n 


X 

n 


F 

— n 


(A  10) 


Since  each  A is  symmetric  and  positive  definite,  it  can  be 

T 

decomposed  uniquely  into  L U , where  L = U is  a lower  triangular 

«n  ,ss  n wn 

matrix  with  positive  diagonal  elements.  Equation  (A10)  is  rewritten  as 


U _ X 


= F 


(All) 


The  solution  for  Xn  can  then  be  obtained  by  first  solving  the  triangular 
system 


for  Z n and  subsequently  solving 


U 


n 


X 


n 


Z 


n 


(A  12) 


(A  13) 


for  X 


n . The  elements  of  /,  are  easily  obtained  from  Eq.  (A  I 2)  in  the  ordet 

»Z  since  lire  first  equation  involves  only  . , the 

2,ii  K,n  l,n' 


-A7- 


< s 


and  so  on.  Ill'  ll  III''  i I < • 1 1 1 « 'ill:; 


second  equation  involves  /. , and  . , 

I , m 2,ii 

ol  X die  similai  ly  obtained  Ituin  llii . (All)  in  the  oidei  X,,  ,X„  , ,...,X, 

n K,ii  K-1,11  l,n 

Those  steps--obtaininq  the  intermediate  solution,  X , Iioim  Eq.  (A  1 2)  and 
obtaining  the  final  solution,  X ^ , from  Eq.  (A13)--are  inferred  tons  forward 
elimination  and  backward  substitution,  respectively. 

The  elements  of  the  symmetric,  pentadinyonal  matrix  A will  be  given 
by 


at- 


The  m.ilrix  1^  ol  the  Lriaiujulai  decomposition  is  given  a: 


dl,n  0 0 


° 1 ,n  d2,n 


fed 
1 1 ,n  2 ,n  3,n 


dK-2,n  0 0 


eK-2,n  dK-l,n  0 


fed 
K-2,n  K-l,n  K,n 


The  elements  of  the  upper  triangular  matrix  are  specified  by  U = 

n 
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at- 


Equaling  the  individual  elements  of  A to  l ho  corre  sponding  elements 

ol  the  product  L U and  solving  lot  the  elements  d,  „ , e,  _ , d„  . 

w n « n 1 , n 1 , n L , n 

l,  . e„  , d„  , etc.  , the  following  relationships  result: 

1 , n 2,n  3 , n 


initialize: 

d 

1 , n 

/ \l/2 

(dl.n) 

Cl,n 

bl,n/‘ll,n 

d2  ,n 

( 2 \ 1/2 
V 2 , n ° 1 » n / 

iterate: 

k 

1,2 K-2 

f, 

k , n 

c /d 

k , n k ,n 

Ck+1, n ' 

( b.  , . - e f 

V k+1 , n k ,n  k 

d,  „ = 

k+2  ,n 

( k+2  , n k+ 1 , n 

1/2 


(A  14) 


Having  determined  the  elements  of  L we  can  solve  explicitly  for 

X in  Eq.  (A10)  by  employing  forward  elimination  to  solve  for  Z in  Eq.  (A12) 
n ~ n 

and  backward  substitution  to  solve  for  X in  Eq.  (A13). 

— n 

The  forward  elimination  proceeds  as  follows: 


initialize:  Z 


l,n 


F,  /d. 

1 , n 1 , n 


2,  n 


F„  - e,  Z 
2 , n 1 , n 1 ,n 


d. 


2 , n 


it<  ratig 


- 3 K 


k,n 


p ( ty  _ r y 

k , n k-  1 , n ‘k-l.n  k-2.n  'k-2.n 
d. 

k , n 


(A  15) 
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I II.  h.i.'kw.iiil  :.uli:;tituli<>n  1 1 1<  ■ 1 1 Iwllow:;: 


initialize:  X„  Z„  /d„ 

K,n  K,n  K,n 


X 


K- 1 , n 


^'K-  I , n ( K-  1 , n AK  ,n 
llK-  1 , n 


iterate: 


K-2 ,K-3 , . . . , 1 


k , n 


'/'k  . n ' k ,n  \tl,n  ' fk  ,n  Xk  + 2,n 
^k , n 


(A  16) 


The  calculation  of  each  element  of  Z of  Z in  the  forward 

K / * * n 

elimination  and  of  each  element  X of  X in  the  backward  substitution 

JS.  f 1 1 I * 

involves  division  by  the  diagonal  element  d . Since  the  system  (A10)  is 

K , 1 1 

to  be  solved  many  times  with  the  same  coefficient  matrix  An  but  with 
different  forcing  vectors,  F , it  would  save  computation  time  if  we  could 
transform  Cq.  (A10)  so  that  the  resulting  diagonal  elements  are  all  unity. 

This  transformation  can  be  accomplished  by  use  of  the  matrices 
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• mil 


■4 


D 


-1 


1/d 


l,n 


1/d 


K,n 


Making  the  transformations 


A = D_1  A D'  1 

an  fa  n » n fa  n 


X = D X 

~ n sj  n ~n 


and 


F = D 1 F 
— n « n ~n 


equation  (A  10)  becomes 


AX  = F 
~ n — n — n 


(A  17) 


(A  18) 


(A  19) 


(A20) 


These  transformations  produce  only  a rescaling  of  the  A ^ , X_n  and  l‘n 
without.  <iffeetin(|  the  character  of  the  coefficient  mail  ices;  i.e.,  A is 

rj  11 

KxK,  symmetric,  positive  definite  and  pentadiagonal . 
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Tin'  1 1 i < 1 1 1 < 1 1 1 1 • 1 1 ( l<  'composi  I Ion  i:.  now  wiilhn 


whore 


and 


Both  L n and  U 


The  elements  of 


a 


b 


c 

where 


from  Lq . (A14)  . 


(A21) 


(A  2 2) 


D 


-1 


(A23) 


have  l's  on  the  main  diagonals  and 


(A24) 


A may  be  written  directly  from  Eq.  (A  17)  as 


k ,n 


a 


k,n/d 


2 

k ,n 


k = 1, . . . ,K 


k,n 


bk,n/dk,n  dk+l,n 


k = 


,K-1 


k,n 


c.  /d 
k , n k 


k + 2 , n 


k = 1 


,K-2 


(A2  5) 
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-AH- 


The  elements  ol  J.  may  be  obtained  directly  from  Cq.  (A 22)  as 


k,n 


k,n 


k ,n 
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-h.,n 


k+i,n 

*k  ,n 

1,  _ 

k+2  ,n 


k = 1,...,K 

k = 1 K-l 

k = 1 K-2 


(A2G) 


Given  the  A matrix,  it  is  necessary  to  compute  and  save  only  the 

d elements  for  the  transformations  (A18)  and  (A19)  and  the  e,  and  f, 
k , n k , n k , n 

elements  for  the  forward  elimination  and  backward  substitution  to  obtain 
X . These  elements  are  obtained  as  follows: 


initialize: 


l,n 


K. 


\ 1/2 


1 ,n 


2 ,n 


bl,n/d!,n 
(a2,n  " el,n) 


1/2 


1 »n 


Cl,n  /d2, n 


iterate: 


k n 


k+ 1 , n 


k t 2 , n 


1,2 K-2 


c /d 

k,n  k,n 


(bk+l,n  " Gk,n  fk,n)/dk+l,n 


/ 2 2 \ 1/2 
V k+2,n  °k+l,n  “ k,n  ) 


k , n 


f,  /d,  „ 
k,n  k i-2,n 


) = o / cl 

k+l,n  k + 1 , n / k+2,n 


(A27) 
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a 


The  d , e ami  1 , elements  may  be  stored  over  the  a,  , b,  and 
k,n  k,n  k,n  k,n  k,n 

c elements,  respectively,  as  they  arc  computed, 
h / * 1 

The  transformed  system  is  solved  for  X by  usiny  forward 
elimination  to  solve  for  Z n of 


Z 

n 


F-n 


(A28)* 


and  backward  substitution  to  solve  for  X n of 


U „ X 


= 5n 


(A2  9) 


The  forward  elimination  proceeds  as  follows: 


initialize: 


‘1 , n 


1,  n 


J2 , n 


* 2 , n ° 1 , n ^ 1 , n 


(A30) 


iterate: 


3 , . . . ,K 


"k,  n 


F.  - e.  . Z.  . - f.  . Z,  . (A3 1) 

k,n  k-l,n  k-l,n  k-2,n  k-2,n 


The  backward  substitution  then  follows: 


initialize:  X^ 

X 


"K , n 


K- 1 , n ?'K-l,n  (K-l,n  \,n 


iti  rale: 


k 


K-2 , K-3 , . . . , 1 


X,  = Z,  - e,  X.  - f . X. 

k,n  k,n  k,n  k+l,n  k,n  k+2,n 


(A3  2) 


From  Eqs . (A  19)  and  (A22)  we  see  that  Z is  equivalent  to  Z 
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could  now  be  used 


The  reverse  transform  of  Eq.  (A  1 0) ; i.e.,  X = D ^ X 
to  obtain  X 

n 

The  preceding  discussion  has  indicated  how  to  solve  Eq.  (A7)  or  its 
equivalent  form,  Eq.  (All)),  by  Gauss  elimination.  Since  this  is  only  one  of 
M equations  of  the  same  form,  which  must  be  solved  for  each  complete  pass 
of  the  relaxation  scheme,  we  will  indicate  how  the  entire  blending  system, 
given  by  Eq.  (5),  is  transformed  as  the  initial  step  in  the  relaxation 
procedure  so  that  line  relaxation  may  be  used. 

The  transformations  (A  1 7 ) , (A18)  and  (A19)  need  to  be  made  only  once 
in  solving  the  blending  system  of  equations  given  by  (5)  or  its  equivalent 
partitioned  forms.  Therefore,  we  will  proceed  by  defining  matrix  D as 


D 

diag  ^ 

2n) 

(A3  3) 

and  its  inverse  D 

-1  by 

Ef1  = 

diag  | 

O'1)  . 

fa  n ) 

(A34) 

Then  (A3)  may  be 

transformed  as 

£_i  ( 

£ + E ) 

D _1  D 

* -1 

P = D R 

(A3  5) 

Let 

D~ ' K 

1 

(All,) 

— * 
p 

P l* 

1 

(A3  7) 

B 

d"1  b 

Dl 

(A3  8) 

« fa 
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r 


and 


A 


D-1  A D 1 


(A  3 9) 


Equation  (A3S)  may  then  be  rewritten 


A P*  + BP 


(A40) 


or  as 


A P*  + B 1 , P*  , + B P* 
s-n  -n  «n-l  --n-l  ^n  — n+1 


= R n = 1 , 2 , . . . , M 


(A41) 


The  solution  by  line  relaxation  for  this  transformed  system  is  defined  by 

(P)(r+1)  - (?„)(r)  - [(P)(rtl/2)  - (?„)(r) 

where  row  ( P*  ) ^ is  the  explicit  solution  to 


(A42) 


L 


n 


p*  \ (r+1/2) 
-n  ' 


(A43) 


by  forward  elimination  and  backward  substitution. 


In  summary,  the  solution  of  (A3)  for  P*  by  line  relaxation  proceeds 
as  follows: 

(I)  Initialize  Lhe  A ami  B mu  trice::  and  the  lotcinq  vector  B . 


( 
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jL 


Decompose1  cuch  A and  save  the  d,  , e,  and  f , 

s.  n k , n k,n  k,n 

elements  u s tn< j liquations  (AX’/).  These  elements  may  be  stored 

over  the  a,  , b,  and  c,  elements,  respectively,  ot  A . 

k ,11  k,n  k ,n  « 11 

Transform  U , I’*  and  B to  k , 1’*  and  B using  Eqs.  (A3t>)  , 

(A 3 7 ) and  (A3U)  . 

Apply  I'.qs . (A1X)  and  (AT 3)  iteratively  until  the  solution  is 
obtained.  During  the  first-half  of  each  pass  the  method  solves 

'k 

explicitly  for  the  1’  associated  witli  the  red  rows, employing 

~k 

the  proper  forcing  function  and  the  previous  pass  values  of  P_n 

associated  with  the  adjacent  black  rows.  The  individual  elements 

of  each  P*  are  then  adjusted  according  to  Equation  (A42)  . During 
~ n 

the  second  half  pass  the  roles  of  the  red  and  black  rows  are 
interchanged . 

Transform  the  solution  P back  to  1’  using  Equation  (A37);  i.e.  , 


aiti:n  i nx  h 


Programming  Details  of  Lino  SOK 


INITIALIZE  A n FOR  n=  1,2,3,...,M 


Diagonal  of  A n 


i = A f B f B +C  +C 

k , n k , n k.n  k + l.n  k.n  k.n-1 


► D.  + D,  *-  E, 

k , n k - 1 , n k , n 


*-  E, 


k-2  , n 


+ F f F + G 

k . n k+2 . n- 1 k . n 


+ G 


k-2 .n-1 


f Rk,n  ' Rk  + 2,n-l  ‘ 4 \ H . i 


+ S.  +•  S + 4 S . 

k , n k-2 . n- 1 k- 1 , n 


+ 16T,  + T 

k . n 


+ T, 


+ T, 


k-1  ,n  k-1 , n<-l  "k+1 , n+1 


+ T 


kU,n 


+ w.  +•  w,  , + w.  , +•  w,  , 

k , n k.n-1  k-i . n k + 1 , n 


for  k odd : k=l,3,5,...,K 


-Bl- 


I 'it  I oil  -i  I i.  l<  |<  >11. 1 1 ol  A 


bk,n  ~Dk,n  2Sk , n 2Sk-l,n-I  4Tk.n  4Tk+l,n  ^k,n  Wk  H , n 


for  k odd:  k - 1,3,...  , K-2 


\,n  = -\,i.n  - 21  Wn-!  - ‘ ’ '"Y.H.n  ' Wk.n-i  ' W 


for  k even : k-  2,4,6, ...,K-1 


Second  off-diagonal  of  An 


’ k , n 


+ 


T 


k * I . n 


d:H,n-l 


W 


k H , n- 1 


for  k odd : k = 1 , 3 , . . . , K-2 


k , n 


-E 


k , n 


k»  1 ,n> 1 


k + l,n 


W 


k ♦- 1 , n 


for  k even : k - 2 , 4 , 6 , . . . , K- 1 


INITIALIZE  I5n  FOR  n 1,2,3 M-l 

Diagonal  of  lin 


o.  = -C.  • T.  . t-  T,  . + W, 

k < n k.n  k-i  ,n  k+1 , n k,  n 


ttk,n  C k , n ]k-l  ,n+l  4 1 k + 1 , nH  Wk,  n 


First:  row  above  diagonal 


\.n 

= 0 

av.,n 

= -rj,  - 2s. 

k , n k , n 

- 4T  - a 

k+1 , n+1 

Second 

row  above  diagonal 

Y„  = 

-G,  + S,  + T 

k,n  k,n 

it 

c 

.V 

-G,  + S.  + T 

k,n  k,n 

First  row  below  diagonal 

^ k ,n 

-B.  , , - 2R  , 

kH,n  k + 1 , r 

- w,  - w. 

u • 1 # n k # n 

Vn  ■ 

0 

41, 


W,  - W, 


2 F. 


4T, 


for  k otJcl:  k=l , 3, . . . ,K 
for  k even : k=2 , 4 , . . . , K- 1 

for  k odd:  k=l,3 K-2 

for  k even:  k=2 , 4 , . . . ,K- 1 

for  k odd:  k=l , 3 , . . . , K-2 
for  k even:  k=2 , 4 , . . . , K-3 

4T 

k , n * 1 

for  I:  odd:  I:  1,3,...,  K-2 
for  k even:  k=2 , 4 , . . . ,K- 1 
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Second  n>w  below  (ii<ni<>ii<>l 

fk,n  " k^2 ,n  ' Kk t 2 , n ' Lk+l,n 

fk,n  ~ ~rki2.n  * ‘V>2,n  f Tk+l,rH  1 


for  k odd:  k=l , 3 , . . . , K-2 
for  k even:  k=2 , 4 , . . . , K-3 
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( 


K.  A p. 

k , n k , n k , n 


R,  b,  + B.  , b.  , . 
k , n k , n k ' 1 , n kH  ,n 


- C.  c,  C,  ,c.  . 

k , n k , n k,n-l  k,n-l 


D,  d,  + D . d - E.  c + E,  e 

k , n k. , n k - 1 , n k - 1 , n k , n k , n k-2,n  k-2,n 


” 1 k . n f k . n ‘ 1 k»-2.n-1  fki2,n-l  Gk , nqk . n ' (‘k-2  .n-lCJk-2  .n-1 


' Rk , n rk t n f Rk+2  , n- 1 rk  <-2  , n- 1 " 2Rk  + l ,n  rk  + l , n 


+ S.  s + S,  _ , s.  0 - 2S  s 

k , n k , n k-2,n-l  k-2,n-l  k-1 ,n  k-1 ,n 


- 4T  t t- T t t- T t +-T  t + T t 

k , n k , n k-1 , n k-1  ,n  k - 1 , n 1 1 k - 1 , n i- 1 k+1 , n^l  k < 1 , n * 1 k + 1 , n k <- 1 , n 

+ W,  w,  + W,  , w,  , - W,  . w - W.  , . w 

k , n k , n k , n- 1 k , n- 1 k- 1 , n k- 1 , n k + 1 , n k t 1 , n 


for  k even:  k = 2,4,6 K-1 
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COM  PUT  i:  n n AND  OFF-DIAGONALS  OF  Ln  TOR  n 1,2,. ...M 


Need  to  save  only  dk  n , , fk  n (write  over  a^  , , c^) 


Initialize: 

dl,n  = 

!=  =»  save 

V l,n 

G1 , n 

bl  .n 

d. 

1 , n 

d2 , n 

Ja„  - e2,  *»  save 

1 2 , n 1 , n 

c, 

1 , n 

e 

1 <—  =»  save 

d2,n 

Iterate:  k = 1,2,. 

. . ,K-2 

ek  + 1 , n 
‘ k ► 2 , n 


°k  + 1 , n 


XlL l 
*k , n 


bk+l,n  °k.n^k,n 
dk+l  ,n 


2 - f2 

ak+2,n  " Ck<-l,n  k,n 


°k* 1 ,n 
dk+2,n 
f,.  _ 


k i-2,n 


save 


save 


save 
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M 


-i  o 


TRANSFORM  R n = D n Rn  TOR  n=  1,2 M 


A 

R. 

R = 
k,n  d. 

k,n 


; k = 1,2, 


I 
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COMPUTE  B n TOR  n = 1,2 M-l 


^.n 


dk,n  dk,n+l 


for  k = 1 , 2 , . . . ,K 


d,  d,  , , 
k,n  k+l,n+l 


for  k = 1 ,2 , . . . ,K-1 


dk,  n dk+2  , n+1 


for  k = 1 ,2 , . . . ,K-2 


6 k , n 


dk+l  ,n  dk,nH 


for  k = 1,2 K-l 


dk+2,n  dk,n+l 


for  k = 1,2 K-2 


TRANSFORM  FIRST  ( 111  ESS  l'(  >R  ( >1)1  > LINES:  .1  1,1 


M 


(PA). 
o l:,n 


<1.  (p  ). 
k , n o 1; , n 


"i  k 1,2,...,  K 


Note:  Tirst  guess  for  even  lines  computed  using  lino  relaxation  (i.e.  , the 
following  steps  with  n even)  with  uu  1 . 


ITERATE  THE  FOLLOWING  STEPS  UNTIL  CONVERGED 

First  half  pass:  n = 1,3,5,  . . . , K 

Second  half  pass:  n = 2 , 4 , 6, . . . ,K-1 


Compute  T n 


k-2,n-l 

Pk-2 ,n-l 

€ k-2 ,n 

P 

pk-2,n+l 

rH 

1 

c 

rH 

1 

l<3a 

— * 

pk-l,n-l 

^ k-1 , n 

“ ic 

Pk- 1 , n+1 

ttk,n-l 

Pk,n-1 

- 

“k,n 

p * 

Fk, n+1 

\,n-l 

Pk+ 1 , n- 1 

5k,n 

— * 

pk+l , n+1 

f k,n-l 

Pk+2 , n- 1 

\.n 

pk+2,nH 

for  k = 1 , 2 , 3 , . . . , K 


Note:  Terms  referencing  indices  outside  of  boundaries  (n  < 1,  n > M, 
k < 1 , k ' K)  are  set  to  zero. 


— B 1 1 — 


;*.ivc  i.n  tin  p;n 


r » 

'<»  P n 


Solve  L n Zn  = Fn 


n 


I nit  ialixe 


Z.  = r. 

1 , n 1 , n 


Z.  = ro  - e,  z. 

2 , n 2,n  1 , n l,n 


Iterate:  k - 3 , 4 , . . . , K 


Z,  = F,  - e,  , Z - f Z. 

k,n  k,n  k-l,n  k-l,n  k-2,n  k-2,n 


Step  2 Solve  U n p*n 


for  Pn 


Initial  ize: 


K-l,n  = Z 

K 1 / n 


* 

eK-l,n  PK,n 


Iterate:  k = K-2  , K-3 2 , 1 


ck,n  PkH,n 


* 

k , n Pk  r-2  , n 


This  becomes  ^P*n^ 


(r+1/2) 
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Correct  a line  of  p*  ^ ; i.c.,  correct  ( p*p 


(^,)<rH,  - (%r  ♦ [(%ri/2)  - (?n  )r;  ■ 


WHEN  CONVERGED 


p*  = p*  /d,  for  k = 1,2,3 K 

*k , n k , n k , n 

and  n = 1,2,3, . . . , M 


IWnmd,  ity  ( '<  >nd  il  n >ns  on  Spread  ing  I’uiamelor  Weights  in  (k.  , n)  Grid  Noluljon 


For  the  following  (k , n)  the  indicated  weights  are  set  : 0. 


1 . First-  1 )if feience  Weights 

k n 


R: 

1 

1,2,3,.. 

. ,M 

first  column  on  left 

1,2,3 K 

M 

top  row 

C: 

2 , 4 , 6 , . . . , K- 1 

M-l 

row  below  top  row 

1,2,3 K 

M 

top  row 

D: 

K 

1,2,3,  .. 

• , M 

first  column  on  right 

1 , 2 , 3 , . . . , K 

M 

top  row 

E: 

K-l  , K 

1,2,3,  .. 

. , M 

rightmost  two  columns 

2,4,6 K-l 

M 

top  row 

T: 

1,2 

1,2,3,  .. 

. , M 

leftmost  two  columns 

2 , 4 , 6 , . . . , K- 1 

M-l 

row  below  top  row 

1 , 2 , 3 , • • • , K 

M 

top  row 

G: 

K-l  ,K 

1,2,3,.. 

. ,M 

rightmost  two  columns 

2, 4, 6,..., K-l 

M-l 

row  below  top  row 

1,2,3 K 

M 

top  row 

II. 

Higher-Order  Spreading 

Parameter  Weights 

R: 

same  as  T 

S: 

same  as  G 

T: 

1,3,5 K 

1 

Dottom  row 

1,2,3 K 

M 

top  row 

1 , K 

1,2,3,.. 

. , M 

left  and  right  columns 

W: 

1 , K 

1,2,3,  . . 

. ,M 

loft  and  right  columns 

1,2,3 K 

M 

top  row 

2 , 4 , 6 , . . . , K- 1 

M-l 

row  below  top  row 
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Spre.tdinq  Parameter  Estimates  in  (k,n)  Notation 


f’stimatc  for 

Estimate  for 

corner  point 

center  point 

Parameter 

(k  odd , m inteqer) 

(k  oven,  m half-inteqer) 

First  Differences: 

, n 

Pk-l,n  Pk,n 

Pk-l,nd  ~ Pk , n 

Ck,n 

Pk  , n *- 1 Pk , n 

same 

C'k , n 

Pkfl , n -p 

k , n 

Pk  t-1 , n+1  ~ Pk,n 

°k , n 

Pk  (-2  , n Pk , n 

same 

^k , n 

Pk-2,nH  Pk,n 

same 

nk , n 

Pkr2,ni-1  Pk,n 

same 

Second  Differences: 

r k , n 

Pk-2  , n+-l  f Pk,n  " 2pk-I  ,n 

Pk-2  , n+I  f Pk,n  " 2pk-l,nU 

st 

k , n 

Pk+2,n+l  Pk  , rr  " 2pk  t- 1 , n 

Pk+2,n+l  f Pk,n  “ 2pkH,nd 

Laplacian: 

1 k , n 

D 4*  D 

k - 1 , n FkM,n 

Pk-  1 , n * Pk+1 , n 

+ n + n 

k-  1 , n- 1 pk  + l , n-1 

f Pk- 1 , n «- 1 * Pk+1 , n+1 

" 4PV-  n 

k , n 

- 4 p 

k,n 

(dross  Difference: 

w 

k,n 

Pk,nH  Pk , n 

Pk,n  + 1 '■  Pk , n 

Pk  - 1 , n Pk  •- 1 , n 

Pk  - 1 , n + 1 Pk  + 1 , n + 1 

-B15- 


Assembly  Module  Coordinator;  (k , n)  of  p and  V . Repor  ts 

i , j i . ) 

Scalar  (k , n) 

p f.  = [(m)R  - - i 

” “ U[()m)r+  0-0r]}t 

First-Difference 

Parameter  (k , n) 

b k = [(j+i)R  - (j-i)J  - 1 

n = U[0“‘)r  + HJ}t 

d k = [((h)t  - (i-i)  J - 1 

n ■ (HHt4-  Hr]}t 

f k . (m)R  - (J-l)R 

n = 'la  [(J+1)r  + Hr"  *]}.,. 

9 k - [Hr  - (j-‘)r  - 2] 

n ’ {i[()H)R  + 0-‘)r-  !Jf.r 

c k = (2i)R  - 1 

p = {HmOd[(2i)Ri2]}t 

o k (21T-l)  MSIGN[MOI){(2j)(<  2},(i-i,r-y 

n - [jt  O.25J  T 


— BIG— 
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